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Multiple Bands — A Key to High - Temperature Superconductivity in 

Iron Arsenides? 
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In the framework of four-band model of superconductivity in iron arsenides proposed by Barzykin and 
Gor'kov we analyze the gap ratios on hole - like and electron - like Fermi - surface cylinders. It is shown that 
experimentally observed (ARPES) gap ratios can be obtained only within rather strict limits on the values 
of pairing coupling constants. The difference of T c values in 1111 and 122 systems is reasonably explained 
by the relative values of partial densities of states. The multiple bands electronic structure of these systems 
leads to a significant enhancement of effective pairing coupling constant determining T c , so that high enough 
T c values can be achieved even for the case of rather small intraband and interband pairing interactions. 



PACS: 74.20.-z, 74.20.Fg, 74.20.Mn, 74.20. Rp 



The discovery of high - temperature superconduc- 
tivity in layered FeAs compounds stimulated active ex- 
perimental and theoretical studies of these new super- 
£h conductors pQ. The main anomaly of these systems is 
their multiple bands nature. Electronic structure in a 
narrow enough energy interval around the Fermi level 
is formed almost only from the d - states of Fe. In fact, 
electronic spectrum of iron arsenides was calculated in 
a number of papers [21 El SI El IS]- The Fermi surface 
consists of several hole - like and electron - like cylin- 
ders and on each of these its "own" superconducting 
gap can be formed. In the energy interval relevant to 
superconductivity electronic spectrum is especially sim- 
ple [HI [9]. It was used by Barzykin and Gor'kov to 
formulate a simple (analytic) model of superconducting 
state of new superconductors [10] . 

Schematically, the simplified electronic spectrum 
and Fermi surfaces of these systems are shown in Fig. 
k> Q] [10]. There are two hole - like Fermi surface cylin- 
^> ,ders surrounding the T point and two electronic pockets 
■around X and Y points in extended Brillouin zone. 

Let Aj be a superconducting order - parameter (gap) 
on the i-th sheet of the Fermi surface. The value of 
A, is determined by self - consistency equation for the 
anomalous Gor'kov Green's function. 

Pairing BCS - like interaction can be represented by 
a matrix: 



C3 



I 

-a 
c 

o 
o 



> 



o 

OS 

o 



V = 



( 


It 


w 


t 


t 


\ 




w 


u' 


t 


t 






t 


t 


X 


M 




V 


t 


t 


f 


A 


/ 



(1) 




2 ' E-mail: 
2 ) E-mail: 



kuchinsk<! 
sadovskiC 



siep.uran.ru 
iep.uran.ru 



Fig. 1. Schematic electronic spectrum and Fermi sur- 
faces of FeAs superconductor in the extended band pic- 
ture. There are two hole - like cylinders around point T, 
while electron - like cylinders are around X (Y) points 

m. 



where matrix elements V l ' J define intraband and in- 
terband pairing coupling constants. For example, A = 
yeX.eX _ yeY,eY determines pairing interactions on the 
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same electronic pocket at point X or Y, /i = y eX - eY 
connects electrons of different pockets at these points, 



V 



hl.hl 



yh2,h2 anc j w _ yh\,hi characterize 



BCS interactions within two hole - like pockets — the 
small one (hi) and the large one (h2), as well as be- 
tween these pockets, while t — V h ' eX — V h ' eY couple 
electrons at points X and T. In Ref. jTU] it was as- 
sumed that u = v! = w. This assumption seems to be 
too strong and below we analyze the general case. 

Superconducting critical temperature T c is deter- 
mined by an effective pairing coupling constant g e ff- 



T r = 



2juj c 



-l/9efl 



(2) 



where oj c is the usual cut - off frequency in Cooper chan- 
nel (assumed to be the same for all types of couplings 
under consideration - a simplification!), while <? e // in 
this model is defined by the solution of the system of 
linearized gap equations: 



where 



9ij 



9eff 
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The matrix of dimensionless coupling constants gij is 
determined by matrix elements of Jl} and partial densi- 
ties of states on different Fermi surface cylinders — Vj 
is density of states per single spin projection on the j-th 
cylinder. 

From symmetry it is clear that v 3 = V4 and the sys- 
tem ([3]) possesses solutions of two types [10] : 

1) solution corresponding to d x i_ y i symmetry, when 
gaps on different pockets at points X and Y differ by 
sign, while gaps on hole - like pockets are just zero: 



Ax = A 2 = 0, A 3 



-A, 



and 



9eff = (p-X)V3. 



(5) 



(6) 



2) solutions corresponding to the so called s pairing 
[3], for which gaps on the cylinders at points X and Y 
are equal to each other: A3 = A4, while gaps on Fermi 
surfaces surrounding the point T are of different sign in 
case of repulsive interaction between electron - like and 
hole - like pockets (t > 0) , and of the same sign for the 
case of t < 0. 

As in this case we have A3 = A4 and v 3 = 1/4, two 
equations in ([3J) just coincide and instead of ([T]), (0| we 



are dealing with 3x3 matrix of coupling constants of 
the following form: 



/ UVi WV2 2tvz 
WV\ u'i>2 2tvz 

\ tV\ tl>2 2\l>3 



(J) 



where A = and ([3]) reduces to the standard problem 
of finding eigenvalues and eigenvectors for the matrix of 
dimensionless couplings gij which has three solu- 
tions, determined by cubic secular equation: 



Det(gij - g e ffSij) = 



(8) 



Physical solution corresponds to a maximal positive 
value of g e ff, which determines the highest value of T c . 

Under the simple assumption of Ref. |TU], when 
u = u' = w, situation simplifies further, as in ([3]) 
only two independent equations remain, so that we have 
2x2 matrix of coupling constants and |(5J) reduces to a 
quadratic equation. Then we easily obtain [TO] : 



Ai = A 2 = kA, A 3 = A 4 = A, 



(9) 



where k 1 = —{g a ff + u(yi + V2)) / '(ti/3) , and maximal 
effective pairing constant is given by: 



2 5e// = ~u{vi + v 2 ) - 2Af 3 + 



(10) 



+y (u(u! + v 2 ) - 2Xv 3 ) 2 + 8t 2 iy 3 ( Vl + v 2 ) 

Possibility of s - pairing in FeAs compounds was first 
noted in Ref. [3]. This kind of solution qualitatively 
agrees with ARPES data of Refs. HU [H EES], except 
the result Ai = A 2 (J9|), which contradicts the estab- 
lished experimental fact — the gap on the small hole - 
like cylinder Ai is approximately twice as large as the 
gap A 2 on the large cylinder. In fact, this contradiction 
is basically due to an unnecessary limitation to the case 
of u = u' = w used in Ref. [TU] . 

The system of linearized gap equations determines 
their ratios on different sheets of the Fermi surface for 
temperatures T — -> T c . In general case, the temperature 
dependence of gaps is determined by the generalized 
BCS equations: 



A, X!-''.' A / 



th 



2T 



(11) 



For T — > these equations take the form: 



A, 



A, 



(12) 
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where we have introduced F(x) = In ( 1+ ^\~ X ~ 

Below we present the results of numerical studies of 
Eqs. ([3]) and (p~2|) for typical values of parameters (cou- 
plings). 

Let us denote the pairing coupling constant on a 
small hole - like cylinder as g = g\\. In the following 
we take g — 0.2, which allows us to remain within the 
limits of weak coupling approximation. 

The ratio of partial densities of states for different 
Fermi surface cylinders in quasi - two - dimensional 
case can be approximated by effective mass ratio on the 
same cylinders. These can be estimated from the data 
for electronic dispersions in symmetric directions in the 
Brillouin zone, obtained in LDA calculations [3 [8], [9]. 
For REOFeAs series (RE=La,Ce,Nd,Pr,Sm...) (1111) 
and for BaFe2As2 (122) from these data we get: 



^ fa 1.18, 
^ ps 1.26, 



^ rs 0.64 
^ ps 0.34 



for 1111 
for 122. 



(13) 



We suppose that pairing interactions on hole - like 
cylinders and between them, as well as on electron - like 
cylinders and between them, arc most probably deter- 
mined by electron - phonon interaction, the relevance of 
which is clearly demonstrated by rather strong isotope 
effect, observed in Ref. [14]. At the same time, inter- 
band pairing interaction between hole - like and electron 
- like cylinders is probably due to antiferromagnetic fluc- 
tuations and is repulsive (t > 0) . It should be noted that 
parameter t from coupling constants matrix ([7]) enters 



Eq. 



determining g e ff, only via t , i.e. independent 



of sign. Thus its sign does not change the value of an 
effective pairing coupling constant and that of T c . Re- 
pulsion between quasiparticles on hole - like and electron 
- like cylinders does not suppress, but actually enhances 
superconductivity leading to the increase of g e ff- Also 
the sign change of t does not change the absolute val- 
ues of gaps on different cylinders, though the repulsion 
between electron - like and hole - like cylinders (t > 0) 
leads to different signs of gaps at these cylinders, while 
for the case of t < both gaps acquire the same sign. 

Despite rather large number of free parameters of 
the model it is not easy to obtain the observable (in 
ARPES experiments of Refs. pH 02 03]) values of 
the ratios |A 2 /Ai| ps 0.5 and |A 3 /Ai| ps 1. In fact 
it requires small enough attraction (or even repulsion, 
v! > 0) on the "large" hole - like cylinder (cf. FigJ2]). In 
the following we assume the ratios of pairing coupling 
constants as w/u = 1, t/u = —1, X/u = 1, which guar- 
antees us the ratio IA3/A1I = 1 for any values of u' and 
arbitrary ratios of partial densities of states at different 
cylinders. Another choice of pairing couplings produc- 



ing IA3/A1I = 1 is also possible, but in general we need 
larger repulsion on "large" hole - like cylinder to get 
IA2/A1I ps 0.5. In Fig. [2] we show the dependences of 
the gap ratios at T = on u'/u, obtained from (fl"2"j) . 
using the partial density of states ratios on different 
cylinders (fl3|) . characteristic for (1111) and (122) sys- 
tems. The gap ratios for T — ► T c differ from the values 
obtained at T = rather insignificantly. 

In Ref. [TS] a two - band model with two hole - like 
cylinders was analyzed, assuming that only interband 
coupling exists, i.e. the coupling constants matrix has 
the form: 



9ij 




(14) 



Under this assumption the gap ratio on hole - like cylin 
ders is given by: 

A 2 _ Pui 
Ai V v 2 



(15) 



so that for characteristic for BaFe^As^ value of ^2/^1 ~ 
1.26 we obtain A1/A2 ps 1.12, which is significantly 
lower than the experimentally observed value of |llj 
A!/A 2 ps2. 

Four - band model somehow similar to that consid- 
ered above was analyzed in Ref. [16j . where tempera- 
ture dependences of gaps (with proper ratios) on dif- 
ferent sheets of the Fermi surface were calculated along 
with the temperature dependence of superfluid electron 
density. However, in this work no analysis was made of 
the important role of multiple bands structure for the 
increase of T c , which we shall discuss shortly. 



-1,0-- ■ 
-1 



— V A i < 122 ) 

A 2 /A, (1111) 

- - - V A , 

g=- UVl =0.2 
w/u=1 , t/u=-1 
(3i+n)/2u=1 



u'/u 



Fig. 2. Dependence of gap ratios on different pockets 
of the Fermi surface on u'/u for g = 0.2, w/u = 1, 
t/u = — 1, X/u = 1 and partial density of states ratios 
given by (Tl3|) . 
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In Fig[3]we show the dependence of an effective pair- 
ing coupling constant and superconducting critical tem- 
perature on u'/u for both classes of FeAs systems (1111 
and 122). It is clearly seen that the effective coupling 
constant g e f / is significantly larger than the pairing con- 
stant g on the small hole - like cylinder. It can be said 
that coupling constants from different cylinders effec- 
tively produce "additive" effect. In fact this can lead to 
high enough values of T c even for relatively small values 
of intraband [2] and interband pairing constants. Actu- 
ally, using this type of estimates we can convince our- 
selves that the critical temperature for superconducting 
transition with d x 2_ y 2 gap symmetry, which is deter- 
mined by an effective pairing constant given by ©, is 
always smaller (for typical values of parameters) than 
the critical temperature for s ± pairing. 




u'/u 



Fig. 3. Dependence of effective pairing coupling con- 
stant on u'/u for g = 0.2, w/u = 1, t/u — —1, X/u = 1 
and partial density of states ratios on different Fermi 
surface pockets given by (|13|l . At the insert — similar 
dependence of the critical temperature. 



To clarify the reasons for the growth of effective pair- 
ing coupling it is helpful to analyze the most simple case, 
when all pairing interactions (both intraband and inter- 
band) in (fT|) are just the same (and equal e.g. to u), and 
all partial densities of states on all four Fermi surface 
pockets are also the same (and equal e.g. to v\). In this 
case we obtain <? e // = 4g = — &uv\, which simply corre- 
sponds to the fact that now the total density of states at 
the Fermi level is four times partial. However, in real sit- 
uation the growth of an effective pairing constant does 
not reduce to this simple summation of partial densities 
of states. In particular, the effective pairing coupling 
may be much larger than the simple sum of intraband 
(diagonal) dimensionless coupling constants, e.g. in case 



of significant interband pairing interaction, which can be 
present in iron arsenides, where the pairing interaction 
between electron - like and hole - like cylinders is most 
probably attributed to antiferromagnetic fluctuations. 

It can be estimated that with the same values of 
interaction constants in ([T]) the critical temperature in 
1111 - type systems is typically larger than in 122 just 
due to the difference of partial densities of states as 
given in (fT5)) (cf. insert in Fig[3]). For example, in case 
of u'/u = (with the values of parameters for 122 - 
system we get the ratio of gaps A2/A1 « 0.6) the cal- 
culated ratio of critical temperatures of 122 and 1111 
systems T c (122)/T c (llll) = 0.67 is very close to the ob- 
served ration of maximal critical temperatures obtained 
for these systems: 38K/55K 0.69. Thus the typical 
difference of T c 's for both classes of new superconduc- 
tors can be attributed to the different values of partial 
densities of states on corresponding Fermi surface cylin- 
ders, despite the fact that total densities of states at the 
Fermi level in these systems are practically the same 
IS]- Of course, the accuracy obtained should not 
be taken too seriously, as in real systems rather strong 
renormalization effects of electronic spectrum (effective 
masses, bandwidths etc.) in comparison with the results 
of LDA calculations are definitely present (and observed 
in ARPES experiments), e.g. due to moderate or prob- 
ably even strong enough Coulomb correlations pQ. The 
main conclusion following from our analysis is the sim- 
ple fact that the value of T c in multiple bands systems 
is determined by the relations between partial densities 
of states on different sheets of the Fermi surface, not by 
the total density of states at the Fermi level as in the 
standard BCS model. 

It should be noted that for the first time (though 
only implicitly) the role of multiple bands structure of 
electronic spectrum as the reason for the increase of su- 
perconducting T c was apparently discussed in relation 
to superconductivity in multivalley doped semiconduc- 
tors jTTJ [TB] . In these works the important role of in- 
terband electron - phonon pairing mechanism was also 
stressed. It was noted that such processes with large 
momentum transfer, leading to reduced screening, may 
be most relevant for the increase of T c . This fact can 
be also important for new superconductors besides the 
abovementioned role of pairing due to spin fluctuations. 

Direct experimental confirmation of the role of mul- 
tiple bands in new superconductors follows from ARPES 
measurements on extremely (hole) overdoped system 
KFe2As2 with T C —3K [19] and similar heavily (electron) 
overdoped BaFe1.7Coo.3As2 [2U], where superconductiv- 
ity is just absent. From these measurements it is clearly 
seen how the disappearance of electronic pockets in the 
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first system and hole - like pockets in the second one 
leads to strong suppression or even the complete disap- 
pearance of T c . 



- (122) 




y — 2a ,a 




---2A 2 /T 




g=-uv,=0.2 




w/u=1, t/u=-1 




(X+n)/2u=1 



-10 12 



u'/u 

Fig. 4. Dependence of ^ ratio on u'/u for 122 - system 
with g — 0.2, w/u — 1, t/u = — 1, X/u — 1 and par- 
tial densities of states ratios on different Fermi surface 
sheets, as given in 1)130 . 

To conclude, on Fig. 0] we show the dependence of 
^ ratio on different sheets of the Fermi surface on the 
ratio of coupling constants u'/u. Here it is important 
to note that the value of this characteristic ratio can 
be significantly different from the standard BCS value 
^ ps 3.5. However, the values shown in Fig. Q] are much 
lower than the ratios observed in ARPES experiments 
[TTl mi IT3"] , where the typical values are 2 ^ 1,3 ps 7.5 
and ^r 2 - « 3.7, which is apparently due to the strong 
coupling effects important in real systems. Our analysis 
was limited to the standard BCS - like weak coupling 
approach. Strong coupling Eliashberg - type analysis 
of multiple bands effects for new superconductors is yet 
to be done. Preliminary results on gap ratios in the 
strong coupling limit for the simple two - band model 
were derived in Ref. [15] . 
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